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Abstract 

Rolls in finite Prandtl number rotating convection with free-slip 
top and bottom boundary conditions are shown to be unstable with 
respect to small angle perturbations for any value of the rotation rate. 
This instability is driven by the horizontal mean flow whose estimation 
requires a special singular perturbation analysis. 

PACS numbers: 47.27.Te, 47.20.Bp, 47.32.-y, 47.20.-Lz 

1 Introduction 

Rayleigh-Benard convection in a plane layer heated from below and rotating 
about a vertical axis, has been the object of special attention motivated 
by both astrophysical and geophysical applications, and by the existence 
of additional instabilities occuring in this system. In the case of free-slip 
top and bottom boundary conditions, Kiippers and Lortz [T] showed, using 
by a perturbation analysis near threshold, that when in an infinite Prandtl 
number fluid, the Taylor number (which measures the rotation rate) exceeds 
the critical value 2285, two-dimensional rolls are unstable with respect to 
perturbations of the form of a similar pattern rotated by an angle close to 
58°. This instability which is also present with no-slip boundaries |2], leads 
in the case of extended systems to the formation of chaotically evolving 
patches of parallel rolls [3], 0], [5], [6], [7]. 

Convection at moderate Prandtl number with no-slip top and bottom 
boundary conditions, was addressed in [8], [2] and [9], and the Kiippers- 
Lortz instability was shown to occur at a critical Taylor number lower than 



1 



in the infinite Prandtl number limit. Free-slip boundaries were considered 
by Swift (cited in [9j) who noted that the usual perturbative calculation of 
the growth rate leads to a divergence in the limit of perturbations quasi- 
parallel to the basic rolls. The present paper is mostly concerned with a 
revisited analysis of this problem, leading to a uniformly valid expression of 
the instability growth rate. We show in particular that for any finite Prandtl 
number and rotation rate, straight parallel rolls are unstable when the angle 
associated to the perturbation is small enough. 

In Section 2, steady convective rolls in a rotating frame are constructed 
perturbatively near threshold. Section 3 is devoted to the computation of 
the instability growth rate for finite angle perturbation, an analysis which, 
at finite Prandtl number, breaks down in the small angle limit. In Sec- 
tion 4, we present a special analysis in the resulting small angle "boundary 
layer", where the interaction of the basic rolls with quasi-parallel pertur- 
bations leads to almost space-independent contributions which become res- 
onant in the zero angle limit. These terms are removed by prescribing a 
quasi-solvability condition the marginal mode of quasi-constant horizontal 
velocity. A uniform expression for the instability growth rate is then de- 
rived and a new "small-angle instability" is obtained. The sensitivity of 
the instability growth rate to the Prandtl and Taylor numbers is analyzed. 
Qualitative features of this instability and its nonlinear development are 
briefly described in Section 5. 

2 Steady convective rolls in a rotating frame 

The Boussinesq equations in a horizontal fluid layer heated from below and 
rotating around a vertical axis z, are written in the non-dimensional form 



where the vertical diffusion time is taken as time unit. We assume a Prandtl 
number > 0.6766, to prevent over-stability [10]. The other parameters 
are the Rayleigh number Ra and the square root r of the Taylor number 
(equal to twice the Rossby number) which, to be speciflc, is taken positive 
(anti-clockwise rotation) . 

Proceeding as in [1], we introduce the operators A = V x (V x .) and T = 
V X ., and express the velocity u = (u, v, u))* in terms of two scalar fields cp 





V.u = 

d 

At} + RaZ.u = u.V?? + — 19 



(2.2) 



(2.3) 
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and -0) in the form u 
where = d^x + d, 
(12.31). we obtain 



yy 



A(0z) + T(-0z) = {dzdx(p+dy'ilj,dzdy(f>-dx'il^,-Ah(f>y 
Applying the operators z.A and z.T on eqs. (j2.ip 
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For free-shp boundary conditions, ?? 
z = ±\. 



dzi^ = in the planes 



A stationary solution of eq. (j2.4p is computed perturbatively near the 
convection threshold by expanding Ra = Rq + eRi + e^i?2 + • • • and X = 
eXi + 6^X2 + 6^X3 + • • • . or, more explicitly, when taking into account the 
boundary conditions satisfied by the individual components. 



(j) = e(f>i cos TTZ + e (j)2 sin 2ttz + • • • 

il) = eipi sin TTZ + £^(-00 + 4^2 cos lirz) + 

1? = ei?! cos TTZ + e^'!?2 sin 27tz + • • ■ . 



(2.5) 
(2.6) 
(2.7) 



Introducing the linear operator L = U +RoG, we get at the successive orders 
of the expansion, 

LXi = (2.8) 
LX2 = -RiGXi + Q{Xi,Xi) (2.9) 
LX3 = -{R1GX2 + R2GX1) + Q(Xi, X2) + Q(X2, Xi). (2.10) 

For a solution in the form of two-dimensional rolls with a critical wavenum- 
ber \ki \ = k, given by the real solution of 



2(^)^ + 3(^)^ 



1 + 



(2.11) 
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the critical Rayleigh number is Rq 
writing, we denote by 



■p |I0]. To simplify the 



Z{a, (3, 7) = (a cos vrz, f5 sin vrz, 7 cos vrz)* 



(2.12) 



vectors corresponding to fundamental modes in the vertical direction and 
obeying the boundary conditions prescribed on X. An element of the null 
space of L is then given by 



ik.x 



v{k) = Z(ci,C2,C3)e 

with ci = 1, C2 = -^, C3 = and kl = + vr^ 

leading order solution reads 

Xi = Av(ki) + c.c, 



(2.13) 



^ , and the 



(2.14) 



where the amplitude A will be determined by a solvability condition arising 
at a higher order. For this purpose, it is convenient to introduce the inner 
product 



{X,X') = Ro / (j)*(l)'dx + Ro / ^p*^p'dx+ / ri^'dx 



(2.15) 



for which the operator L is self-adjoint. 
Using the notation 

Q(X„ X,) + QiX„X,) = (1 + 5.,)(Qg, qSJ, Qg)*, 



(2.16) 



we have in eq. (|2.9 
Q(Xi,Xi) = 



2i2o7rf^sin 27rz\A\'' 



(2.17) 



The solvability condition for eq. (j2.9|) (obtained by taking the inner product 
of this equation with v{ki)), requires Ri = 0. Defining the operator P = 
+ r'^dzz - Ro^h, eq. I^M reduces to 



(2) 



)(3) 



(2.18) 



The right hand side vanishing identically, we get Ah(j)2 = 0, since elements 
of the null space of P, already included in (pi, are not needed in (/)2. 
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For the two other components of X2, one easily checks that -02 = 0, 
^0 = m^e^ik^-^ + c.c. and i?2 = ©2 with ^-q = ^^A^ and 62 = ^I^P- 
This enable us to compute 

/° - - ^ 

Q{Xi,X2)+Q{X2,Xi)=\ ^ ^0 (^*e**^i-^-3^e3**^i-^) sin 7rz + c.c. . 

\ 27rfc^02^e*'^i-^ cos vrz cos 27rz + c.c. / 

The solvability condition of eq. (j2.1Up then reduces to R2 = r2\A\'^ or 
equivalently, e^l^P = , with r2 = ohi(Rok^ ~ "^bt")) which completes 

the computation of the roll amplitude in terms of the distance to threshold. 



3 The Kiippers-Lortz instability 

We assume that the steady rolls of wavevector ki computed in Section 2 
are subject to a perturbation X in the form of rolls with an infinitesimal 
amplitude and a wavevector k2 making with ki an angle 9 that it is enough 
to consider in the range ] — f , f ]• We assume for the sake of simplicity that 
the wavenumbers \ki \ and \k2\ are critical. When real, the growth rate a of 
this perturbation is given by 

{U + RaG)X = Q{X, X) + Q{X, X) + aVX. (3.1) 

In order to compute a = eo\ + e^cT2 + • • • perturbatively near threshold, we 
also expand X = X\^ eX2 + e^^s + • ■ ■ or, for the individual components, 

(j) = (pi cos TTZ + e(j)2 sin 27rz + • • • (3-2) 
ip = ipi sin TTZ + e(0o + ^2 cos 27rz) + • • • (3-3) 
"!? = i?i cos TTZ + ei?2 sin 2ttz + • • • . (3-4) 

Equation (13. ip leads to 

LXi = (3.5) 

LX2 = QiXuXi) + Q{XuXi) + aiVXi (3.6) 
LX3 = Q{X2,Xi) + Q{Xi,X2) + Q{Xi,X2) + Q{X2,Xi) 

+aiVX2 + (T2VX1 - R2GX1 . (3.7) 

Using a notation similar to (j2.16p . we define 

Q{X„X,) + Q{X„X,) = Q^j = (Qg,Qg,Qg)*. (3.8) 
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Writing the solution of eq. (|3.5p in the form Xi = Bv{k2) + c.c, where B is 
an arbitrary constant, we have in the right hand side of eq. ()3.6p . 



j^{5^AB*e'''~-^ + S^ABe'''^-^) + c.c. 
^ h{5+AB*(^^~-^ + 5^ABe^^^-^)siii2TTz + c.c. J 

where we have introduced the wavevectors 

fc^ = fei zb A:2, 

and defined the numerical constants 



3 ^2 . 



Furthermore, the coefficient 6^, given by A^e 



TTT' 



1,2 



Rok'^TT 



2kl 



(3.9) 



(3.10) 



(3.11) 



, read 



5+ = -4fe2 cos2 - , 67 = -4^ sin2 - . 
h 2 h 2 



(3.12) 



Since {v{k2),Q{Xi,X{) + Q{Xi,Xi)) = 0, while {v{k2),VXi) / 0, the 
solvability condition for eq. (j3.6|) implies cJi = 0. Straightforward algebra 
then leads to 



A#2 = Qf| + i2oA/,<A2. 



Solving in the form 



(3.13) 
(3.14) 
(3.15) 

(3.16) 



4>2 


= ^+ABe'''^-^ 


+ ^^AB*e'''~ 


+ c.c. 


(3.17) 


^2 


= ^AB(^^^-^ 


+ ^^AB*e'''~ 


'■^ + c.c. 


(3.18) 


4^0 


= ^+ABe^^^-^ 


+ §0 ^5*e*^" 


'■^ + c.c. 


(3.19) 


^2 


= Q+ABe^^*-^ 


+ 62 ^5*e^^" 


'■'^ + c.c., 


(3.20) 



we get 



6 



''■4. _ _ -^'■"h^ _ ^ 

©2^ = i?0^(7&f'^+ + JS) + JsJ^ ©2" = i?0^(i^5- + JS) + Jsl^ . 

The coefficient 6^ andp^ defined by the condition Ae'^^^'^T{27rz) = 5^e^'^^-^T{2TTz) 
and Pe^^^'^T{2TTz) = p^e^^^'^T(2TTz), (where the function T stands for sine 
or cosine), are given by 



<5+ = - 


-[47r2 + 4A:^ cos^-], 


(3.21) 


(5- = - 


-[47r2 +4A;2 sin^^], 


(3.22) 


P^ = - 


9 
-[47r2 + Ak^ cos^ -f - Att^t^ + 4k'^Ro cos^ -, 


(3.23) 


P~ = - 


-[Att^ + Ak'^ sin^ - Att^t^ + 4/c2i?o sin^ ^, 


(3.24) 


•? ^ and 


sin^ ^ contributions result from the action of the hori- 



zontal Laplacian on g*''^-^ and e*'^ respectively. 

An important observation is that the contribution 'i'^ AB*e^^ to ^2 
(which disappears at infinite Prandtl number) diverges in the limit 6 — > 
0, where it can be viewed as associated to a "mean flow" generated by 
the rotation. This term is specific to free-slip boundary conditions and 
has no equivalent when rigid boundaries are considered. The divergence 
originates from the fact that in eqs. (I3.13p - (l3.16p . the dynamics of the 
mean flow is slaved to that of the leading convective mode. This "adiabatic 
approximation" is valid at finite 9 but breaks down in an "angular boundary 
layer" near ^ = 0, where time derivatives become relevant. Postponing to 
Section 4 the analysis of this layer, we derive here the solvability condition 
of eq. p.7p for finite 9, by writing 

{v{k2),GXi) = ^ciCsk^B* (3.25) 

(vih), VXi) = ^{RoPr-'cjk^k/ - RoPr-^cie + 4)3* (3.26) 

{v{k2).{Q{X2,Xi)+Q{Xi,X2)) = -^c\clk^\A\''B\ (3.27) 
Furthermore 

g(Xi,X2) + Q(X2,Xi) = Z{q%q%qfl)\A\^B*e-'^^-^ + c.c. + F,, (3.28) 
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where F3 refers to non resonant terms proportional to sinSyrz or cosSvr^;. 
We also have 

gj^ = — sm6'[(-^+ + 2^o)(ci7r2-ci/c2-2ci7r2cos6l-2c27rsm6l) 

+(^2 - )(ci7r^ - ciA:^ + 2ci7r^ cos6' + 2c27rsin6') 
+2l>^A;p(c2(l - cos 61) + vrci sin 61) - 2^^A;p(c2(l + cos6l) - vrci sin^)] 

(3.29) 

= - — [2§o(ci7r(-l + 2cos6'-cos26') + C2(sin6'-sin26')) 
J^'^ 2iv 

+2^'J"(ci7r(-l - 2 cos 61 - cos 26*) - C2(sin + sin 261)) 

sin6'(-2ci7rsin6' + C2(l + 2 cos 6*)) 

sin6l(-2ci7rsin6l - 02(1 - 2 cos 6*)) 

(C27r(l - cos 29) - vr^ci (2 sin 9 - sin 261)) 
+$2 (C2vr(l - cos 261) + vr^ci (2 sin 6^ + sin 26^))] (3.30) 



/t2 



^-^2 2^^^ (ci7r(l + COS0) + C2 sin6') + 62 (ci7r(l - cos0) - C2 sin6')] 

+ ycsl^-^ + 21-+ - ^-2 - 21-0 ] sin 0. (3.31) 
It follows that 

{v{k2),Q{Xi,X2) + Q{X2,Xi)) = i(iiocigg + i2oC2(?3 + C3(?g)|A|'i?*, 

(3.32) 

and finally 



2 _ ^r2CiC^k'-{^cicik^ + R^Ciql!^+RoC2q\l+c^ql!^) ^ 

~ ' RoPr-^clk%^ - RoPr-^4k^ + 4 ' 

(3.33) 



where e^|Ap can be expressed as ^^^^a^4 {Ra — Ro)- 



Since in the limit ^ 0, ^'g diverges like sin ^ |, the quantity gj'^ with 
/ = 1, 2, 3, scales like sin ~ sin~^ |, and the growth rate behaves like 

^^2^21^12 g2 

indicating a breakdown of the above asymptotics at finite Prandtl numbers, 
in the case of small angle perturbations. 
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Pushing the ^-expansion at the next order, (as needed in Section 4), we 
write 

^ ~ [-^ + S^(2e + (3.35) 



with 



and 



2/c|Pr- sm I 



^ r2_ 

"^o(l + i(l-Sf?)) 



(3.36) 
(3.37) 



the latter coefficient collecting contributions originating from . 

The divergence shown in eq. (j3.34p was noted in [9] . It indicates that the 
above analysis should be viewed as an outer expansion, and that a different 
scaling is required for small 9. In the following, the growth rate given by eq. 
()3.33p will thus be denoted (t outer- 

4 The small-angle instability 

The small angle divergence of the stream function -00 e sin~^ | and of the 
growth rate Couter ^ sin~^ |, indicates that new scalings in e are expected 
in an angular boundary layer near ^ = 0. Denoting by e° the thickness of 
this layer, by the amplitude of ^'o and by e"^ the magnitude of the growth 
rate in this layer, the matching of the "outer" and "inner" regions requires 
(3 = 1 — 2a and 7 = 2 — a. Since, in the inner region, the time derivative in 
the mean flow equation (whose presence will remove the divergence) becomes 
comparable to the viscous term when 7 = 2a, we get a = |, /3 = — ^ and 

7 = 1- 

Furthermore, when expanding eq. (j2.4p inside the boundary layer, the 
parameter e appears not only through the horizontal Fourier modes of Xi 
whose amplitudes scale like entire powers of e, but also through the angular 
dependence of the operators involved in this equation. We are thus led to 
expand 



and 



4 5 n 7 8 

a = eai + e3(j4 + e^as + e (T2 + e^a? + eaas + ■ ■ ■ (4.1) 
33 33 



X = e ^Y_i +Y0 + X1 + e^Yi + e^Y2 + eX2 + etj^i + etXs 

3 33 33 
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(4.2) 



where terms of the form = ^ ^B* Ae^'^ '^(0, 1,0)* + c.c, are introduced 
to cancel almost resonant contributions resulting from the interaction of the 
basic rolls with quasi-parallel perturbations. As seen later, in the boundary 
layer, a can be complex. 

Substituting (14. 2p and (14. ip in eq. (j3.ip and concentrating on perturba- 
tions such that the angle 6 between the wavevectors ki and k2 is of order 

2 

e 3 , we obtain the following hierarchy. 
• At order e'^, 

LXi = 0, 



leading to 

• At order e, 
( 

LXo = 



Xi = v{k2)B + c.c. 







+ (TiBVv{k2) + c.c. 



\ -4j3k'^AB*e''' •^sin27rz 
The solvability condition reads 

ecTi = 0, 

and the solution is given by 



/ 



Xo 





^^ABe^''^-^ 



+ c.c. 



with ^+ = and 02 



V 62 AB*e''' ■^sin27rz 

fc2 



33-. 



At order e 3 . 



elLX, = ehme^*^,\A\^Be''^^--^Z{-^ 

3 3 



+e3a4VBv{k2) + c.c. + MTZ, 



where AfTZ collects non-resonant terms. The solvability condition is 



£3 0-4 = e/c^ sin 

3 3 



(4.3) 
(4.4) 

(4.5) 



(4.6) 



(4.7) 



(4.8) 



(4.9) 
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• At order e 3 , 

3 

+eia5VBv{k2) + c.c. + MTZ. (4.10) 
The solvability of this equation requires 

ehi=ePsme\A\^^*Q. (4.11) 

• At order e^, 

+e3sin0§lie^'=2-^|^|25Z(-2— C27rsin0 , — ci7r(sin - - sin — ) , 0) 

3 ^^T" -^^V ^ ^ 



e^R2GXi + e^a2VBv{k2) + c.c. + A/'Te, (4.12) 



where 



iQ%Q%Qfl) = (0,4:-ci7r^+,A:^ci7re2-) (4.13) 

denotes the limit as ^ ^ of the vector {q^^ , , Q^^^ ) from which the 

contributions coming from ^'q have been removed. The solvability condition 
reads 

= -i]\A\'^e'^ + 2ehk'^ sine sin -\A\'^^Zl + el k'^ sine\A\'^^l* , (4.14) 

23 3 

where 

r2CiC3fc2 - (icfcife^ + Ihciq\'l + RoC2q?l + o^qf-) 

f] = f b£ i2± b±. (4 15) 

RoPr~^4k'^k/ - RoP-^clk'^ + ^ ' ' 

identifies with the expression given in eq. ( I3.36P . The coefficient ^ is given 
by eq. ^M- 

Combining the solvability conditions (|4.6p . (I4.9p . (I4.11|) and (j4.14p . we 

are led to express the growth rate 

(dinner = eSfJi + etfJs + e^(T2, (4-16) 

3 3 



k^ 
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in terms of the "mean flow' 



^ = e 1 + §n + e^§i 



(4.17) 



3 3 



in the form 




sme\Af{l + 2^sin-)^' - e^r]\A\'^ 



(4.18) 



where subdominant corrections have been neglected. 

In order to estimate the mean flow ^, we push the e-expansion of eq. 
()2.4p at higher orders, where the beating of the perturbation with the basic 



independent and thus resonant in the small 9 limit. Consequently, uniform 
boundedness of the solutions requires, in addition to the usual solvability 
conditions, the prescription of "quasi-solvability conditions" aimed to elim- 
inate terms which are strictly resonant only for 9 = 0. This approach is 
similar to that used by Ablowitz and Benney [11] when dealing with the 
small-amplitude divergence of the Whitham modulation analysis for nonlin- 
ear dispersive waves (see also [E]). These authors modify the (algebraic) 
dispersion relation by means of additional corrective terms determined by 
a constraint which becomes an actual solvability condition in the small am- 
plitude limit, thus transforming the algebraic dispersion relation arising in 
Whitham's theory, into a partial differential equation for the wave ampli- 
tude. In the context of rotating convection, we include contributions Yn 
proportional to ^'-i^i-^'i)-^ in the perturbation expansion, which are deter- 
mined by cancelling them with the terms displaying the same functional 
dependence and originating from the beating of the basic rolls with quasi- 
parallel perturbations. Like in the small-amplitude limit of nonlinear waves, 
this condition becomes an actual solvability in the limit 9^0. In both 
instances, the singularity is prevented by removing slaving conditions: that 
of the amplitude with respect to the phase in the case of waves, or that of 
the mean flow with respect to the convective modes in the present problem 
(compare eqs. ()3.14p and (|4.25p below). 

7 

• At order e s , 



solution produces contributions of the form e*^'^^ ^'^^'^ which become space- 



/ 6(^,7 e^^B*e*^ ■^sin27rz \ 



3 " -3 



V J 



ei[g^4] + e ^J^AB*e'''~-^ + C.C. 

\ js6j^ABe^^'^-^ sm2Trz + c.c. j 



+ 



12 



+ 



36^ 



esm9^'^e"'^-''\A\^BZ{-2—C2TTsm9 , — ci7r(sm - - sin — ) ,0) + 

3 -'T" I'^y 3 3 

(4.19) 

where [Q i] denotes the leading order of Q{Xi,Xa) + Q{Xa,Xi). Although 

X± contains terms proportional to e^*''^'^, the resulting contributions of the 

3 

form e^^ in Q 3 are preceded by a factor proportional to sin^ f and thus 

1'3 ^ 

not included in [Q 4]. The quasi-solvability condition thus reads 

^ 1 "o 



0. 



(4.20) 



• At order es 



e^LXs + 

3 



V 





3 3 3 I -,i 





e2 A fc^^-o (5e- 3 e**^ + c.c.) sin e cos' 







v 







+ 



1,4: 1,4: Q oa 

€3 sine§*ie*'^2-^|^p5Z(-2^C27rsin0 , — Ci7r(sin - - sin — ) ,0). 

3 Pr Pr 2 2 



(4.21) 



The quasi-solvability condition is 

-2 T 



- p-h-^a,fi% + 5J^'% - P-'eia,_6^{e-^>*_^) = (4.22) 

where, as previously, [Qi ] does not contribute. 
• At order e^, 



e^LXi = e3[Q^ 3 + Q3 J + - €^RiGXi - e^i^sGXa + £3(^2^X2 + a^VXi) 
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-(0 ,e3 5^^^tAS*e~*^ ■^ + c.c. ,0)* 



+(0 ,P~'6f^e-^{a2^*_^+a4m+a,%)AB*e~'^ + c.c. ,0)*, (4.23) 

3 33 3 

with the quasi-solvabihty condition 

6fie-^^^*)-P-^e^a2Sj;{e'H*_^)-P-hlai6^{eH\)-P-hla,6j;% = 0. 

3 3333 

(4.24) 

Combining eqs. (I4.20p . (I4.22p and ()4.24p . we get, up to subdominant 
contributions, 

- P-V„^ + ^h'^ = ^Pr^h (4.25) 

which together with eq. (I4.18p . constitute a closed system. Solving the 
resulting quadratic equation for the growth rate, we obtain two solutions 

4nner = ^("eM^P " 4A;2p, sin^ ^) ± 



(e%|Ap - Ak'^Pr sin^ -f - Ae^k^ s\.u6\A\^ i2{l + 2 sin 



(4.26) 



where rj, ^ and j2 are defined by eqs. ()3.36p . ()3.37p and ()3.1ip . This expres- 



sion covers several regimes 
fi) For 6* ~ €3 , 



~ -{-Ak^PrSin^ - ) ± - 



{ik'^Pr sin^ -)2 _ 4e2^2 sin6'm2j2 



(4.27) 

In this range, (jf^^^^ > for ^ > (0 < 0) if r > (resp. r < 0) for any 
finite value of the Prandtl number (still assuming Py > 0.6766) and of the 
Taylor number. 

(ii) When 6* > ei, 

^tnner ~ ^match = [-^l + ^^{^C + -^)]eM\^ (4-28) 

iK^rr sm I 

and matches the limit of cTcmter as — s- 0. Similarly, 

^rnne.~-4A;'^.sin2^ (4.29) 

is negative and becomes of order unity outside the boundary layer. 
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(iii) For 6* ~ 



^tnner ~ ^("^'r^l^p) ± 1 [(e^T^I^H^ - Ae'^ sme\A\''h\ ' . (4.30) 

and for 6 = 0, cr+„g^ vanishes, while cr^nneri^) = -e^??|^P- 

We thus obtain a uniform representation for g] — ^ , ^] , of the instabihty 
growth rate near the convection threshold, of the form 

(^'^ = Crtiner + Pouter - <7 match- (4-31) 

where the various terms arising in the right-hand-side of eq. ()4.3ip are given 
by eq. p.33p . (|4.26p and ()4.28p . The influence of various parameters like 
the Prandtl number and the rotation rate on the stength of the instability, 
is illustrated in the following figures. 

Figure 1 shows the variation of the eigenvalues with the angle 6 of 
the perturbation for = 2, e = 0.1 and r = 10. For anti-clockwise rotation 
and finite Prandtl number, the growth rate is positive for small enough 
positive angles 9. There is also a range of negative angles, where there are 
two complex conjugate eigenvalues, with negative real parts. The dashed line 
represents the outer solution a outer which diverges in the limit — > 0. The 
other eigenvalue (j~ which, as e ^ 0, becomes marginal in a neighborhood of 
= 0, is of order unity outside the angular boundary layer. It thus cannot 
be computed perturbatively for 6 order unity but, being always negative or 
complex with a negative real part, it cannot lead to an instability. 

Figure 2 displays the growth rate for r = 38, e = 0.1 and various 
values of the Prandtl number for positive angles. We observe that both 
the range of unstable angles and the maximal growth rate decrease when 
the Prandtl number is increased. At P^. = 10, the small angle instability 
and the Kiippers-Lortz instability (around 9 = 50°) can be separated, in 
contrast with the case of smaller Prandtl numbers (e.g. Pr = 5) where all 
the angles < ^ < 64° are unstable. For this rotation rate, only the small 
angle instability survives at Prandtl number P^ = 15. It becomes hardly 
visible at Pr = 50. Indeed, as the Prandtl number goes to infinity, the 
negative eigenvalue a~ has a limit, while the outer expansion o-^te^ extends 
towards = where it asymptotically reaches the value (T~(0), the inner 
range reducing to the vertical axis. 

Figure 3 shows the variation of the instability growth rate with the ro- 
tation rate r, for P^ = 15 and e = 0.1. For r = 10, only the small angle 
instability is present. The Kiippers-Lortz instability (again localized around 
9 = 58°) arises for r ~ 40 and is strongly amplified as r is increased. 
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Figure 4a displays for e = 0.1, the critical value of the rotation rate r 
for the onset of the Kiippers-Lortz instability, as a function of the Prandtl 
number, as long as the latter is large enough for the two instabilities to be 
separated. Figure 4b shows the most unstable angle (in degrees) for the 
Kiippers-Lortz instability, versus the Prandtl number, for a rotation rate 
corresponding to the onset of the instability. 



16 



0.5 i 




Figure 1: Instability growth rates and cr^ or their real part r when 
complex conjugate (full line), together with the diverging "outer solution" 
(dashed line), versus the perturbation angle 6 (in degrees), for P = 10, 
r = 10 and e = 0.1. 
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Figure 3: Growth rate , versus the perturbation angle ^ > 0, for P = 15, 
e = 0.1 and different values r = 10, 40, 60, 100 of the rotation rate r. 




Figure 4: Kiippers-Lortz instability boundary in the (P^-, ''")-plane (a), and 
angle associated to the unstable perturbation at the critical Taylor number 
versus the Prandtl number (b). 

5 Nature of the instability and nonlinear develop- 
ments 

We showed in Section 4 that in a rotating horizontal fluid layer with mod- 
erate Prandtl number, limited by top and bottom free-slip boundaries, con- 
vective rolls are linearly unstable with respect to perturbations in the form 
of rolls making a small angle with that of the basic pattern. This instabil- 
ity occurs even when the rotation rate is too low for the existence of the 
Kiippers-Lortz instability. It is related to the divergence of the growth rate 
()3.34p which, at finite Prandtl number, occurs when the direction of the 
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wavevector of the perturbation, approaches that of the basic rolls. We here 
chosed the associated wavenumbers to be critical, but the effect survives 
whatever their values. 

The above instability was obtained in an infinite domain. Its persistence 
with lateral boundaries requires the presence of a large number of rolls, 
and thus a convective cell with a large aspect ratio The mesh size in 

Fourier space scaling like /i, the minimum angle 9 between two wavevectors, 
behaves like /i-^/^. Since for the small-angle instability a ~ e^/^ and 9 ~ e^/'^, 
it follows that a ^ fi, a, growth rate intermediate between that (of order 
unity) of a pure amplitude instability and a phase instability, for which a 
scales like /i^. 

As shown in [18] and |19j in the absence of rotation, that of parallel rolls 
may also be unstable (for wavenumbers larger than critical) to a skewed- 
varicose instability whose growth rate also varies like the inverse aspect ratio 
/i, a scaling resulting from the strong magnitude of the mean flow in the case 
of free-slip boundary conditions. This instability is however not captured by 
the present formalism since the roll distortions involved in this instability 
cannot be represented within the class of perturbations (superposition of 
two families of straight rolls), we have considered. 

In order to investigate the relation between the small-angle and the 
skewed-varicose instabilities, and to analyze their nonlinear developments, a 
system of equations in the spirit of the Swift-Hohenberg equation, but cou- 
pling the leading vertical mode to the mean flow, was derived by a systematic 
perturbation expansion near threshold [13]. This system which preserves 
the rotational invariance of the problem, generalizes equations obtained by 
Manneville [H] at finite Prandtl number in the absence of rotation. In a sim- 
plified version where the non-local couplings are suppressed and only a few 
representative nonlinear terms are retained, it is also consistent with models 
used in |15j and [16] for rotating convection at infinite Prandtl number. A 
similar model was considered in |17j . 

As discussed in [13], the phase equation derived in the context of the 
generalized Swift Hohenberg equations, shows that the skewed varicose in- 
stability occurring without rotation near onset, becomes asymmetric with 
respect to the angle of the phase perturbation, in the presence of rotation. 
This model also shows that both the asymmetric skewed varicose and the 
small angle instabilities lead, by means of reconnection to a progressive ro- 
tation of the convective rolls in the direction of the external rotation, an 
effect of the mean flow which develops shear layers. 

We are thus, led to conclude that the small-angle divergence of the 
Kiippers-Lortz instability growth rate pointed out in [9], results from the 
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presence of a small-angle instability, which can be identified as an asym- 
metric skewed- varicose instability. In contrast with the (symmetric) skewed 
varicose instability arising without rotation, the asymetric one developing 
in presence of rotation exists whatever value of the basic roll wavenumber. 
Both instabilities are produced by the mean flow and deaseapear in the limit 
of infinite Prandtl numbers. 
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